function [theta, phi, R] = TotalAntennaResponse() 
phi = linspace(pi-pi/200, -pi, 200);
theta = linspace(pi,0, 200);
phip = zeros(200*200,1);
thetap = zeros(200*200,1);
detNames = {'INDIGO'};
R = zeros(200*200,1);
 for k = 1:length(detNames);
     Det = LoadDetectorData(detNames{k});
    for i = 1:200;
        for j = 1:200;
            thetap(200*(i-1)+j) = theta(i);
            phip(200*(i-1)+j) = phi(j);
                [Fp , Fc] = ComputeAntennaResponse(phi(j), theta(i), 0, Det.d);
                R(200*(i-1)+j) =  R(200*(i-1)+j) + sqrt(Fp.^2 + Fc.^2);
        end
    end 
    
        figname = ['fig/' detNames{k} '.png']
        load coast;
        figure(8);clf;
        
       
    projectionType = 'mollweid';
    mapstruct = defaultm(projectionType);
    mapstruct.angleunits = 'degrees';
    mapstruct.origin = [0 0 0];
    mapstruct.falseeasting = 0;
    mapstruct.falsenorthing = 0;
    mapstruct.scalefactor = 1;
   
    
    xproject([thetap,phip], sqrt(R), projectionType);
    colorbar;
    caxis([0.0, 2.0]);
    xlabel('\phi / rad', 'FontSize', 14);
    ylabel('\theta / rad', 'FontSize', 14);
    title(detNames{k}, 'FontSize', 16);
    set(gca, 'FontSize', 12);
    box on;
    print('-dpng', 'detNames{k}''.png')
    
    hold on;
    % load the worldmap coastlines:
    load coast;
    [xc, yc] = mfwdtran(mapstruct, lat, long);
    plot(xc, yc, 'k-', 'LineWidth', 0.2);
    
    % Plot the locations of all the detectors on top of this
       [x, y] = mfwdtran(mapstruct, Det.phi, Det.lambda);
       plot( x, y, 'ro', 'MarkerFaceColor', 'r', 'MarkerSize', 8);
       hold off
 end
    
   
       
%     figname = ['fig/LHO+LLO+Virgo+INDIGO.png'];
%        
%     projectionType = 'mollweid';
%     mapstruct = defaultm(projectionType);
%     mapstruct.angleunits = 'degrees';
%     mapstruct.origin = [0 0 0];
%     mapstruct.falseeasting = 0;
%     mapstruct.falsenorthing = 0;
%     mapstruct.scalefactor = 1;
%    
%     figure(4); clf; 
%     xproject([thetap,phip], sqrt(R), projectionType);
%     colorbar;
%     caxis([0.0, 2.0]);
%     xlabel('\phi / rad', 'FontSize', 14);
%     ylabel('\theta / rad', 'FontSize', 14);
%     title('LHO+LLO+Virgo', 'FontSize', 16);
%     set(gca, 'FontSize', 12);
%     box on;
%     print('-dpng', 'LHO+LLO+Virgo+INDIGO.png')
%     
%     hold on;
%     % load the worldmap coastlines:
%     load coast;
%     [xc, yc] = mfwdtran(mapstruct, lat, long);
%     plot(xc, yc, 'k-', 'LineWidth', 0.2);
%     
%     % Plot the locations of all the detectors on top of this
%     for j=1:4
%        det = LoadDetectorData( detNames{j} );
%        [x, y] = mfwdtran(mapstruct, det.phi, det.lambda);
%        plot( x, y, 'ro', 'MarkerFaceColor', 'r', 'MarkerSize', 8);
%     end
%     hold off;
%     
    
         
        
